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Abstract 

Using the Fock-Schwinger proper time method, we calculate the induced 
Chern-Simons term arising from the Lorentz- and CPT-violating sector of 
quantum electrodynamics with a b^tp^^ip term. Our result to all orders 
in b coincides with a recent linear-in-6 calculation by Chaichian et al. [ |hcp 



th/0010129l v2]. The coincidence was pointed out by Chung [Phys. Lett. 
B461 (1999) 138] and Perez- Victoria [Phys. Rev. Lett. 83 (1999) 2518] 
in the standard Feynman diagram calculation with the nonperturbative-in-6 
propagator. 
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The Lorentz- and CPT-violating Chern-Simons modification to the Maxwell theory was 
first proposed a decade ago An important feature of the Chern-Simons term || is 
that Lagrange density is not gauge invariant, but the action and equations of motion are 
gauge invariant. In their Lorentz violating extension of the standard model, Colladay and 
Kostelecky posed the question whether such a term is induced when the Lorentz- and CPT- 
violating term ^frfgip (b^ is a constant four vector) is added to the conventional Lagrangian 
of QED i- 

Recently, several calculations have been carried out to determine the radiatively induced 
Chern-Simons term from the Lorentz- and CPT-violating fermion sector P|HT5|]: 

£ = - e jjL - 75J 6 - m)V> ■ (1) 

Jackiw and Kostelecky clarified that the induced Chern-Simons term is not uniquely deter- 
mined ||. It depends on whether one uses a nonperturbative definition or a perturbative 
definition of the theory defined by Eq. In a nonperturbative formalism, radiative cor- 
rections induce a definite nonzero Chern-Simons term, while when a perturbative formalism 
is used, radiative corrections are finite but undetermined. The regularization scheme one 
chooses to adopt can generate further ambiguity in both nonperturbative and perturbative 
formalisms. 

In the standard Feynman diagram calculation with the nonperturbative-in-6 propagator, 
Chung U and Perez- Victoria |7| demonstrated that the result to all orders in b coincides with 
the previous linear- in-6 calculation by Chung and Oh Q as well as Jackiw and Kostelecky || . 
More recently, nonstandard approaches were employed also in the calculation of the induced 
Chern-Simons term. Chan |T(J and Chaichian et al. |12], used the covariant derivative 



expansion Jl7 and the Fock-Schwinger proper time method Hl8,19fl, respectively, to obtain 



the induced Chern-Simons term. The common feature in these two nonstandard approaches 
is to develop a series of local effective Lagrangian in powers of II M = id^ — eA^, rather than 
in powers of id^ and separately. 

The purpose of this work is to calculate, using the Fock-Schwinger proper time method, 
the induced Chern-Simons term arising from the Lorentz- and CPT-violating sector of QED 
with a ipfyy^ip term keeping the full b dependence in order to see whether the coincidence of 
all-order-in-6 result with linear-in-6 calculation takes place in this nonstandard approach. 

The effective action, r eff , of the theory defined by Eq. (|l|) is given by 

r eff = -i\n(JV'ipVil;e i f d4xC ^ = -iTr ln($ - - -y 5 p - m) , (2) 

where the trace "Tr" is taken over both spinor indices (of any combinations of 4 x 4 Dirac 
matrices) and space-time coordinates. 

Let us decompose the trace in Eq. (Q) in the following manner, as was done in Eq. (3) 
of Ref. pfl: 



1 In the nonperturbative definition, we use the 6-exact propagator S(p) = for the standard 

perturbation calculation. (See the work of Colladay and Kostelecky [16| for detailed analysis of 
the extended Dirac theory with this propagator.) In the perturbative definition, the propagator 
has the 6-independent form S(p) = and —75^ is considered as an interaction vertex. 
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Tr ln(z^ — eJji. — 75/$ — m) = Tr ln(i^ — e^L — m) 



■ 1 ^ Tr f — - 1 75^ ) • (3) 



Then, the effective action can be written down as follows: 



- _ r (o) , r (i) 

eff — 1 e ff eff 1 



where 

,(0) 



T e g = —iTi ln(i^ — ejfi. — m) , 

i$ + e4 + zjj+^l 1 ^) ' (4) 

The induced Chern-Simons term is contained in r^. By neglecting the dependence of b 
in the denominator of the trace in , the authors of Ref . [0] obtained the undetermined 
induced Chern-Simons actionP] 

TqI = ^ J d 4 x e^^bpAvFxp . (c being an arbitrary contant) 

If one keeps the full dependence of b in the calculation of for extracting the induced 
Chern-Simons term, one expects that the induced Chern-Simons action would take the 
following form: 



r 



(i) _ ce 2 L J b 2 



cs- 4^2 



1 + / 



J d A xe^ Xp b^A v F Xp , (5) 



where f(b 2 /m 2 ) is some function of a single argument b 2 /m 2 . In what follows, we are to 
show by explicit calculation that this function f(b 2 /m 2 ), indeed, vanishes. 

Now, let us introduce a (fermionic) Green's function G(x, x') as the inverse of the operator 
— i$+eA'(x) + zj 5 p+m, i.e., let us assume that G(x, x') satisfies the following inhomogeneous 
differential equation 

[-z# + e4(x) + z>y 5 p + m]G(x, x) = 5 4 (x - x) . (6) 
Then, T^g of Eq. (|j) can be written as 

T$ = J d A xC§ = -iJ tfx j\lztT[G{x,x') lb f>] x ^ x . (7) 

The trace "tr" in this equation is taken over the spinor indices only and the limit x' — > x is 
taken by averaging the forms obtained by letting x' approach from the future and from the 
past [11J. Further, if one introduces a (bosonic) Green's function A as follows: 



This undeterminicity arises from an intrinsic ambiguity in the limit linL E _ + o x p x v /x . (See Eq. (36) 
below.) This limit has directional dependence in a strict mathematical sense as emphasized in 
Ref. ill. 
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G(x, x) = — eA\x) — zjsp + m]A(x, x) , (8) 

or, equivalent ly, 

G(x, x') = [-id' p - eA^x')] A(x, x')^ + A(x, x') \-z lb p + m] , (9) 
then, one finds that Eq. (|j) becomes 

HA(x,x') = 5\x-x') , (10) 

where Ti is defined as 

U = -n 2 + m 2 + z 2 b 2 + ^uF^ + 2iza^b vlb , (11) 

with II M = id p — eA p . Borrowing the quantum mechanical matrix element notation, the 
Green's function A(x',x") can be expressed as follows: 

POO 

A(x',x") = (x'\H- l \x") =i ds (x'\e- iHs \x") . (12) 

Jo 

In the s integration, a convergence factor —it (e > 0) in 7i is understood. The idea of the 
Fock-Schwinger proper time method is to consider 7i in Eqs. ( |TUD and (|TT| ) as a Hamiltonian 
that governs the evolution in "time" s of a hypothetical quantum mechanical system. Then 
the integrand in Eq. (|12"D becomes a transition amplitude (x'(s) \x"(0)). Thus, the solution 
A(x', x") can be written, in terms of this transition amplitude, as follows: 

/•oo 

A{x',x")=i ds(x'(s)\x"(0)) . (13) 
Jo 

Now, let us evaluate the transition amplitude (x 1 \s)\x" (0)) . It satisfies the evolution 
equation 

id s (x'(s)\x"(0)) = (x'(s)\n(x(s),U(s))\x"(0)) , (14) 

with the boundary condition 

(x'(s)\x"(0)).^ = 5\x'-x"). (15) 

In order to obtain the Hamiltonian in the Heisenberg representation, 7i(x(s), Il(s)), one first 
has to know the evolutions of x and II. Using the commutation relations 

[z'*, rr] = -ig^ , [if, ir] = -ieF^ , 

one can obtain the equations of motion for the operators x(s) and Il(s): 
Xf, = i[H, Xp] = 2I1 M - 2iza fiu b v j 5 , 

n M = i[H, n M ] = 2eF pv W + iedTFp, + \d v F pv aP l/ - 2izeF iip a pv b ulh . (16) 
In the constant F pv approximation, these equations are solved, in matrix notation, as follows: 
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n(s) = e 2eFs (n(0) - izab-ys) + izab-y 5 , 

x(s) = x(0) + (eF)' 1 ^ 8 - l)(n(0) -izabj 5 ) 

From these two relations, one can find 

n(s) = ieFe eFs [sinh(eFs)]- 1 [x(s) - x(0)] +^675 , 

n(0) = ^eFe- eFs [sinh(eFs)]- 1 [x(s) - x(0)] + iz<r&7 B 
n 2 (s) = [x(s) - x(Q)]K[x(s) - x 



+ize[x{s) - x(0)]e- eFs [sinh(eFs)]" 1 Fa67 5 - 3z 2 b 2 



where 



(17) 



(18) 



[IV 



K = ^e 2 F 2 [smh(eFs)}- 2 

Further, using the commutation relation 

[x,(s),x u (0)]=i(eF)-\e 2eFs - 

one can fincf] 

x(0)Kx(s) = x(s)Kx(0) - ^TR[eFcoth(eFs)] . 

From Eqs. (|TT|), (p~8"l), and (|19|), one finally obtains 

W(x(s),n(s)) = + 2x(s)ifx(0) - x(0)A"x(0) 

-^TR[eFcoth(eFs)] - |TR(aF) + m 2 - 2^ 2 6 2 

Therefore, Eq. fll4"|) becomes 

z\9 s (x'( S )|x"(0)) = 

whose solution is 



(a;' - x")AT(x' - x") - -TR[eFcoth(eFs)] 



-TR(aF) + m 2 - 2z 2 b 2 



(x\s)\x'\0)) , 



(x'(s)\x"(0)} = C(x',x")s~ 2 e- L ^exp 



- ~{x' - x")eFcoth{eFs){x' - x") 



le. 



+ ^TR(aF)s - i(m 2 - 2z 2 b 2 )s 



(19) 



(20) 



3 The trace denoted by "TR" extends over Lorentz indices only: TR(AB) = A^ U B U/1 . In particular, 
note that the expression TR(cri ? ) = a^ u F UIM is still a matrix valued quantity because each is a 
4x4 matrix. 
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where 



L(s) = -TRln[(eFs)- 1 sinh(eFs)] . 
From the definition of the operator II and the first two relations in Eq. flT8|), one can find 

^eF[coth(eFs) + 1](V - x") + izabj 5 ) (x' \s)\x" '(0)) , 



[id 1 -eA(x')}(x'(s)\x"(0)) = 
[-id" - eA(x")](x'(s)\x"(0)) 
from which, in conjunction with Eq. (0), the differential equations for C(x', x") are obtained: 



-eF[coth(eFs) - 1] {x - x") + izabj 5 ) (x'(s) \x"(0)) , (21) 



id' - eA(x') - -F{x' - x") - izab^ 



-id" - eA(x") + -F(x - x") - izab-y 5 



C {x j x 
C{x', x h 



0, 
. 



Two forms of the solution for C(x',x") are obtained: 



C(x',x") =d(x")exp 



z{x — x")ab^5 — ie j dx fJ '(A(x) + —F{x — x" 
J x" ^ 2 



C(x r , x") = C 2 (x') exp z{x - x")ab l5 + ie f dx^(A(x) + -F( 

L Jx' ^ 2 



(22) 
(23) 



Since both integrals in Eqs. ( p2|) and ( p3|) are independent of the integration path due to 
the Stokes' theorem, one may choose the integration path to be a straight line connecting 
x' and x", i.e., x = x" + t(x! — x"). Then, it is readily shown that 

f dx»(F(x - x")) = C dt(x' - x")F(x' - x") = , 

Jx" V ' i 1 JO 

£ dx»(F(x - x'))^ = J° dt(x' - x")F(x' - x")(l - t) = . 

With these vanishing integrals in mind, the comparison of two expressions for C(x',x") in 
Eqs. (|22|) and (|23| ) leads us to the following conclusion: C\(x") = C^x') = C (constant). 
The constant C is finally determined by Eq. (|T^) as 



(4tt) 2 ' 

Therefore, the complete form of the transition amplitude is given as 



(x'(s)\x"(0)) 



(47T) 

ie. 



; $(x',x")s- 2 e- L ^exp 



-~(x' - x")eFcoth(eFs)(x' - x") 



+ -TR(aF)s - i(m 2 - 2z 2 b 2 )s 
2 



(24) 



where 
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z(x' — x")ahy 5 — ie dx^A^{x) 



$(x', x") = exp 



Now, let us compute tr[G(x, x^j^p] in Eq. (0). From Eqs. 



(25) 

and (p~3|) , we obtain 



tr[G(x,x')j 5 p] = tr 



(i^-^K) / ds(x(s)\UJs)\x'(0)) 



+izb tr 



tfa<a;(s)|a/(0)> 



tr 



ilJf-l^K) / ds(x(s)|n„(0)|a;'(0)> 



+226 tr 



ds(x(s)\x'(0)) 



Averaging these two equivalent expressions gives 



tr[G(x,x') l5 p] 



-tr 



75/ ^(x( s )|[n M ( s )-n M (o)]|x'(o)) 



6,, 

tr 

2 



750" 



+izb 2 tr 



ds(x( s ) l + n u (o)]|a/(o)> 

ds(x(s)|z'(0)) 



(26) 



From the first two relations in Eq. (|18D , we have 

(x(s)|n M (s) - n M (0)|x'(0)) = [eF(x - x')Ux( S )|x'(0)> , 

(ar(s)|n M (s) + n M (0)|x'(0)) = [eFcoth(eFs)(x - x') + ^0-6 75 ] M (:r(s)|a/(O)) . 

Using Eqs. p^ ) and (p?D, we finally rewrite Eq. fl26| ) as follows: 

1 



(27) 



tr[G(x,a;')75^ 



(4tt)< 



» (is 
2" 



$ (X,X') / ^ g-ifm'-a^Vjs e -L(s) ^ifx-i'JeFcothfeFs)^-!') 



S z 



7s ^(z-z')^ e f TK(aF) s 



ie 



+ ^F fia [coth(eFs)} a P(x - x') p b u ti 



ixv e z(x-x')ab e fTR(aF)s 



-b 2 tr 



e z(x-x')ab e fTR(aF)i 



(28) 



where $0(2^)^') is the function Q(x,x') in Eq. ( |2"5] ) at z = 0. 

The vacuum current vector (J^x)) associated with 6 is obtained from by variation 
of ^(2): 



6T 



(i) 

eff 



(29) 



Thus, from Eqs. (0) and 
positive real axis, 



8A»{x) ' 

we have, after deforming the integration path of is to the 
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le 



;(X - x')^ (x,x') [ dZ [ % 

Jo Jo s z 



(4tt) 2 

„-(m 2 -2z 2 b 2 )s _-^(s) „i(x-x')eFcot(eFs)(a;-a;') 



X 



ieF p ,(a; - x') v Vtr 



7s 



e 2(a!-a:')<T676 e f TR(aF) S 



+ -F pct [cot(eFs)] Q/3 (x - x')A tr 



T5Cr P y e z(x-x')ub l5 e f TR(ctF) S 



where 



1 



-TRln[(eFs)- 1 sin(eFs)] . 
2 



Using the eigenvalue technique Hl9| , e ^ is determined as 

(esfg 



Im cosh(esX) 



where 



X = ^2(F + ig), T= - V" = — -(E - B< 



£ = -F^F^ = -E 



with F MZ/ = ^t pvXp F\ p . From the matrix decomposition formulas 

e z(x- X >)aby 5 =ai + a2 ( x _ x ') a b l5 ; 

e f tr( ct f)s = Ci + C2TR ( orF ) + , C375 + iC4 TR(aF) 75 

where 



oo ^ 2 JV 



ai= L(m 



[{x-x'fb 2 -({x-x')-b) 2 ] N , 



oo ~2AM-1 



° 2 £(2iV+l)! 



[(ar - rr') 2 & - ((x - x') ■ b)T , 



ci = Re cosh(esX) , c 2 = Re [sinh(esX)/(2X)] 
c 3 = Im cosh(esX) , c 4 = Im [sinh(esX)/(2Jf )] 



the trace quantities in the integrand of Eq. ( p0|) are evaluated as follows: 

tr [ e *(*-*'>&7 6 e fTR(<xF )s j = 4fliCi _ g ia2C2 ( x _ x ^P b _ Si a2C ^ x _ ^JTft ; 

tv c-' "'" ' <>" 1 = 4iaac 3 - 8a 2 c 2 (x - x')Fb + 8a 2 c 4 (x - x')Fb 

-8i ai (c 2 F pv + c A F pu ) 



tr 



75(7^ e ^ a '-^)^75 gf TR(«tF) S 



+ 4a 2 ( Cl [(x - z')^ - (x - z'pl + ic 3 (x - if) a bp^ Ppv ) 
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+ 8ia 2 (c 2 [b p ((x - x')Ff - b v ((x - x')F) p 

+ {x- x'yibFf -{x- x') p {bF) v \ 

+ Aic 4 [e pua(3 {x - x') a {bF) p - e puaf3 b a ((x - x')F) p 

- b p ({x - x')F) v + b v ((x - x')F) p 

+ {x- x r ) p (bFf - (x - x'YibFY}) . (32) 

Now substituting Eqs. fl31"|) and (|32|) into Eq. fl3~0|) and discarding terms that do not 
contribute to the Chern-simons current, we obtain 



(J^x)} = — (x - x)^ (x,x') I dz I ds 



1 



o Jo Im cosh(esX) 

xe -(m 2 -2z 2 b 2 )s e l(x-x')eFcoth(eFs)(x-x') 

x (a lC2 eF pa [cot(eFs)] a/3 (a; - x)pb u F pp + a 2 c 2 zb 2 (x - x') p F pv b v 



(33) 



Thus, the current (J M (x)) of Eq. (|33|) is cast, in the leading order approximation in the 
coupling constant e, to the following form:[] 



,2 ,1 



In 2 jo 



dz 



(x-x')^x-x') p %(x,x')F pu h 



x ds[ — + 



a% a 2 zb 2 



-(m 2 -2z 2 b 2 )s+(x-x') 2 /(4s) 



(34) 



In the light of 6-perturbation theory, the insertions of ^75 are revealed best via the expan- 
sion of the factor e 2z2h2s [hi Eq. fl34|) 1 in powers of its exponent. Then using an integration 
formula 



with K n (x) being the modified Bessel function of order n, the Chern-Simons current (J^x)) 
in Eq. fl34|) is evaluated as follows: 



,2 ,1 



2tt 2 



dz 







( / >"» rp ^ ry ^ 



■6^ 



00 2V 



X £(2iV)!l 1 + 2iV + 



(x — x') 2 

/) 2 \ 

" T — 2 [-m 2 (x-x') 2 ] 1/2 ) 
1 m z / 



\l '\2r2 ( X x )a(x X W/ t\2lxn~B 

x< (x — x 7 ^ -(x — x ) b tr 

I [x — x') 1 

00 2n , 1.2 x n 

E ( [-^(^ - x') 2 ] ( " +1)/2 Xn-i [-m 2 (x - x') 2 ] 1/2 



n=0 



(35) 



This approximation matches with the conventional Feynman diagram calculation of the one-loop 
vacuum polarization tensor. 
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Using the following formulas 



lim $o(x, x ) = 1 , 

x'—*x 

lim [-m 2 (x - x') 2 ] 1/2 KM-m 2 (x - x') 2 ] 1/2 ) = 1 



\im[-m 2 (x - x') 2 )( n+1 V 2 K n J[-m 2 (x - x') 2 ] 1 / 2 ) = 0, n=l,2,3,--- 

lim — — x x ^ = C g ) ( c being an arbitrary contant) (36) 

x'-*x [x — X'Y 

it is not difficult to see that each higher order term separately vanishes. Thus, we finally 
arrive at 

ce 2 ~ n 



W*)> = -^ F ^y dz 



and equivalently 



cs 



ce 2 

471-2 



which means that the all-order-in-6 calculation does not alter the linear-in-6 result and thus 
the functiom f(b 2 /m 2 ) in Eq. (|5|) vanishes. This complete our calculation. 

We understand the identity of the lowest-order calculation with the all-order cal- 
culation by following an argument of Coleman and Glashow [pOj] . Since in the expansion of 
the 6-dependent vacuum polarization amplitude in powers of b (6-perturbation theory), all 
except the first order are free of linear divergences. Hence there is no ambiguity in evaluating 
higher order graphs. Momentarily let each of the two photons carry different momenta, say 
Pi and P2 (this means that the chiral inserions of ^75 carry non-zero momentum); Applying 
the Ward identity to this vacuum polarization amplitude, one finds that the amplitude is 
0(pi) and 0{j>2) at the same time, i.e., it is O{pipo); now go to equal momenta, p\ = P2 = P, 
and observe that the amplitude must be 0(p 2 ). The Chern-Simons term one is seeking is 
0(p); hence all these higher-order graphs do not contribute. 

In summary, we have used the the Fock-Schwinger proper time method to calculate the 
induced Chern-Simons term arising from the Lorentz- and CPT-violating sector of quantum 
electrodynamics with a b^ip^^^^ip term. Although there is an ambiguity in the calculation as 
mentioned in footnote 2, our all-order-in-6 result coincides, independently of this ambiguity, 
with the linear-in-6 result. 
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